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Abstract 



The effect of angular momentum conservation in microcanonical tlier- 
modynamics is considered. This is relevant in gravitating systems, where 
angular momentum is conserved and the collapsing nature of the forces 
makes the microcanonical ensemble the proper statistical description of 
the physical processes. The microcanonical distribution function with 
non-vanishing angular momentum is obtained as a function of the coor- 
dinates of the particles. As an example, a simple model of gravitating 
particles, introduced by Thirring long ago, is worked out. The phase di- 
agram contains three phases: for low values of the angular momentum L 
the system behaves as the original model, showing a complete collapse at 
low energies and an entropy with a convex intruder. For intermediate val- 
ues of L the collapse at low energies is not complete and the entropy still 
has a convex intruder. For large L there is neither collapse nor anomalies 
in the thermodynamical quantities. A short discussion of the extension of 
these results to more realistic situations is exposed. 



1 Introduction. 



Conventional thermodynamics applies to systems whose forces saturate and 
therefore macroscopic parts have negligible interactions. If this is the case, 
the macroscopically conserved quantities are extensive, the thermodynamic po- 
tentials are homogeneous functions of degree one of these quantities and large 
systems can be studied in the thermodynamical limit. These properties hold if 
the so called stability condition is verified. For a system of N classical particles 
interacting through a two body potential, (t>{r), the stability condition states 
that there must exist a positive constant Eq such that, for any configuration 
{ri, . . . , r^}, the following inequality is obeyed 



$(r-i,...,rAr) = -^^(|r,-r,|) > -NEo . (1) 

These cases constitute very specific situations. For example, (||) holds for po- 
tentials repulsive enough at short distances, like the Lennard- Jones potential, 
but it is clearly violated for gravitating systems. When the stability condition 
does not hold, the system undergoes a phase transition from a high energy gas 
phase to a collapsing phase at low energies. This low energy regime is not a 
proper thermodynamical phase, since it is not homogeneous, but, nevertheless, 
we shall use the name collapsing phase throughout this paper. 

Microcanonical thermodynamics has been recognized as the statistical de- 
scription of systems which suffer fragmentation and clustering [Q, since there 
seems to be problems to describe spatial inhomogeneities within the canonical 
formalism. These phenomena, fragmentation and clustering, appear in many 
different branches of physics, from nuclear physics and atomic clusters [Q 
to astrophysics § . The task of microcanonical thermodynamics is to compute 
the entropy of a given system as a function of the macroscopically conserved 
quantities. It is intuitively clear that the entropy of a system which undergoes 
a phase transition associated to spatial inhomogeneities like clustering or col- 
lapsing will depend crucially on the value of its total angular momentum, if 
conserved, since the spatial distribution will do. (See ref. for a similar dis- 
cussion in a different context). Angular momentum is of particular importance 
in astrophysics: the formation of a single star or of a binary system, the birth of 
a solar system around a star or the merger of galaxies in clusters is, to a large 
extent, determined by the value of the angular momentum. 

It is well know that, as a consequence of the virial theorem, gravitating 
systems have negative specific heat [R, |1. In fact, it seems that the common 
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feature of all these collapsing systems is the appearance of an interval of energies 
where the microcanonical specific heat is negative Since the canonical 

specific heat is always positive, the microcanonical ensemble gives the proper 
statistical description of these systems. 

A system of classical gravitating particles interacting via the newtonian po- 
tential has an infinite entropy, due to both short and long distance singularities 
. It is clear that the short distance singularity of this potential is not physi- 
cal, since new physics (like quantum mechanics) appears at small scales. Hence, 
there should exist a natural short distance cut-off. The long distance singularity 
has a different nature. In principle, one would say that gravity should be able 
to bind the system. A little caveat, however, easily convinces oneself that, at 
least from the statistical point of view, the system prefers to evaporate than to 
remain bound. One needs a box (the long distance cut-off) to keep the system 
confined, but no such a box appears in nature. We have to consider it as an 
artifact making the statistical description of the system possible, which will be 
sensible if the evaporation rate is small. The box breaks translational invari- 
ance, and therefore the momentum is not conserved. Since we are interested in 
keeping angular momentum conserved, we must deal with spherical boxes, in 
order to maintain the rotational symmetry exact. 

The paper is organized as follows: in section 2 we compute the microcanon- 
ical distribution with conserved angular momentum, by integrating out the mo- 
menta in the phase space. The formula obtained gives a microcanonical weight 
suitable for Monte Carlo simulations, and allows also for analytical approaches 
relying in mean field methods. In sec. 3 we briefly derive the mean field equa- 
tions for systems which conserve angular momentum. Sec. 4 is devoted to the 
discussion of a simple model, introduced by Thirring which mimics the main 
features of gravitating systems with great success. The model is extended to 
take into account the conservation of angular momentum. The article ends with 
a summary of the conclusions, in sec 5. 



2 The microcanonical distribution. 

Consider a system of N classical particles whose interactions are given by 
a general potential energy, <I>, depending only on the position of the particles. 
The hamiltonian reads: 
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If the system is isolated and the potential energy is translationally and ro- 
tationally invariant, the energy, momentum and angular momentum will be 
conserved, and, as a consistency condition, the center of mass will move with 
constant velocity. Without any loss of generality, we can take the total mo- 
mentum zero and then the center of mass is fixed, for instance, at the origin of 
coordinates. Assuming some kind of ergodicity, the microcanonical distribution 
will give the statistical description of the system. The entropy S is given by the 
Boltzmann formula: 



SiE,L,N) = log WiE,L,N) , (3) 

where E is the total energy and L the modulo of the angular momentum, since, 
by symmetry, the entropy must be independent of the direction of the vector 
L; W is the volume of the phase space shell defined by the orbits with given 
E and L. Notice that we took the Boltzmann constant equal to unity, and the 
entropy dimensionless. Hence, the temperature is measured in units of energy. 

To avoid too cumbersome expressions, let us consider the case in which only 
the angular momentum is conserved. If also the linear momentum is conserved, 
the derivation of the microcanonical distribution is similar and the result will 
be reported at the end of the section. The volume of the relevant cell of the 
phase space can be computed starting from its definition: 



W{E,L,N) = ^ I n^^^ S{E-n)S('HL-Y.^.^P.). (4) 

Integrating out the p variables, we will obtain a non-singular microcanonical 
probability distribution depending only on the spatial configuration {ri, . . . , rjv}- 
In this way, we shall get an expression suitable for microcanonical Monte Carlo 
simulations. 

To perform easily the p integration in (^) , let us define 



/jv 2 
2—1 i i 

where we used the notation E ^ E — ^{ri, . . . ,rN). It is convenient to take the 
Laplace transform of 11 respect to E, 
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poo 

n(s,L,7V,M) = / dEe-'^Ii{E,L,N,{r}) , (6) 
Jo 

with Res > 0. Introducing the foUowing representation for the remaining Dirac 
delta: 



K^) - ^e-, (7) 

we obtain 



t\{s,L,N,{r}) = 



(8) 



The integral in p is now gaussian and can be readily performed. What shall 
remain is again a gaussian integral in u;; there is no difhculty in evaluating it. 
After some algebra it is found: 

n(s, L, iV, {r}) = C (det I)-^'^ \[^r!^jyj^ . (9) 

where I~^L = X]a/3=i I^p ' ^^"^ matrix / is the inertial tensor respect 
to the origin, 



N 

4=1 



3jY — 3 3/2 

with a,/? = 1,2,3 labeling the coordinates, and C = (27r)~2— Hi "^i is a 
constant. 

The inverse Laplace transform of (^) can be found in any book of integral 
transform tables It is: 
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if o_ < £; < 6 



(11) 



Therefore 



3N-5 



Il{E,L,N,{r}) - ^;p|^(det/)-i/2 (^e - ^L'^I'^L ] " (12) 



ii E > ^ L^I ^L, and it vanishes otherwise. 

After the p integration the volume of the phase space is given by: 




W{E,L,N) I [\\^] (i? - Il^I-'L $ 



(13) 



where C = C/[iV!r((3iV — 3)/2)]. The thermodynamical quantities can be 
computed numerically from the last integral by, for instance, using a suitable 
Monte Carlo algorithm. Eq. (|l^) can be also used to derive microcanonical 
equations for mean field approximations. 

If the linear momentum is also conserved, the same expressions hold, chang- 
ing the exponent in eq. (O) by (3A^ — 8)/2, the constant C by 



C = (2.)^(;^m,.)-/^n-''' (14) 



and the argument of the T function appearing in C by (37V — 6)/2. In this case, 
the inertial tensor / refers to the center of mass, which can be taken at the 



origin, and then eq. (10) holds. Notice that, if the number of particles is large. 



the momentum conservation gives negligible differencies. 

If the particles have some internal spin, then the intrinsic inertial moment of 
each particle must be added to the orbital inertial moment, and the exponent 
in (^) must be properly modified to take into account the number of intrinsic 
degrees of freedom. The derivation of a microcanonical weight for these more 
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general cases is identical to that outlined in this section and presents no addi- 
tional difhculty. As a last remark, we stress that the results of this section, i.e. 
eq. (^3|), apply to stable as well as to unstable systems. 



3 Mean Field theory. 

If the number of particles N is very large, one can expect that the force that 
a particle undergoes will be more sensitive to the mean particle distribution of 
the system than to the fluctuations around it. Then one can think that the 
thermodynamics of the system depends only on the mean particle density. In 
this section, we shall give a derivation of the mean field equations using this 
hypothesis (cf. rcf. fl^). For reasons which will become clear in the following, 
the results of this section apply only to unstable systems. 

Let us start with the expression given by eq. (|l3|). To perform the integral 
in r, we divide the integration region in A cells of volume a^. The integral over 
(firi becomes a sum over all possible cells. We have then N of such sums. For 
each configuration {r} we can define a function n{i), with i = 1, . . . , A, which 
counts the number of particles in the i-th cell. It is not difficult to see that the 
initial sum over the position of the particles can be reorganized as a sum over 
the number of particles in each cell, as follows: 



where ^4 is a constant which contains C, the elementary phase space volume 2iTh 
and the volume of the cell, a^. The combinatorial factor arises since a certain 
amount (= this factor) of configurations of particle positions give the same 
cell occupation distribution {n{i)}. Introducing the particle density p{ri) = 
n{i)/{Na^), where denotes the position of the center of the i-th cell, we can 
write (p^) as the following functional integral: 




3JV-5 




(15) 




[dp]S 
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exp<^ N 



- J d\p{r) [logp(r) - 1] + ^ log L^I-^L - $ 



(16) 



where we have substituted the factorials by its asymptotic form since we consider 
N large, and ignored some irrelevant constants. We have also neglected the term 
VdetT, which is of order 1/N . Now, we have to clarify how the argument of 
the logarithm scales with A''. For a given discretization, it is clear that the 
potential energy <I> scales as iV^, irrespective of whether the potential is stable 
or unstable. With N fixed, in the continuum limit, a — > 0, the potential energy 
scales with N if the potential is stable, and a proper thermodynamical limit 
can be taken. For unstable potentials, however, the potential energy scales 
with N"^ still in the continuum limit. The usual thermodynamical limit does 
not exist for this kind of systems. However, we can take N ^ oo and send 
the masses and the coupling constants in the potential energy to zero, keeping 
Nm = M and N'^<j){r,r') constant. In this case, the energy E is not scaled; 
the thermodynamical functions depend on the total energy, instead of on the 
energy per particle. However, as can be seen from (|6|), the entropy scales with 
N. This is the scaling to be considered in this paper, and applies naturally to 
unstable systems. 



The functions / and $ in (|l6|) depend on p{ri) through a suitable discretiza- 
tion. For very large N and very small a, they can be expressed as: 



I^fi[p] = M J d^r p{r) {r^5^p-r„rp) 

= j d^rd^r' (j){r,r')p{r)p{r') . (17) 



The functional integral (^^ is defined by the previous discretization. There 
are two limits to be considered in (^6|), a ^ and N oo. This is the 
rigorous order of the limits. For unstable systems, mean field theory consists 
in interchanging both limits, taking first — > oo. In this case the functional 
integral (|l6|) is saturated by the maximum of the exponent in the integrand, 
all the fluctuations around the maximizing density being suppressed by powers 
of 1/N. Ignoring all the irrelevant constants, we can define the entropy per 
particle as 



S = - j d?rp{r) [logp(r) - 1] + ^ log (^E - \l^I-^L - <i>^ , (18) 
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with / and $ given by (p^. The physical distribution p is that which maximizes 
S under the constraint Jp(r) = 1. 



The distribution which maximizes (18) is one for which a principal inertial 



axis goes along the angular momentum direction, for, if not, one can define a new 
distribution by rotating the given distribution until the principal inertial axis 
with larger inertial moment coincides with the angular momentum direction. 
The potential energy <i> and the pure cntropical term J p{[ogp — 1) do not 
change, but the rotational energy has been thus lowered, and hence S increases. 
Taking the angular momentum along the z axis, we can write 



5 = - y" d?rp{r) [logp(r) - 1] + | log (^E 



$ ) . (19) 



2L 



33 



The maximum can be obtained by taking the functional derivative respect to p, 
which gives the following integral equation: 



p{r) = exp<{ 



a^' + y') - / dV0(r,r')p(r') 



(20) 



where p, is the Lagrange multiplier for the constraint J p = 1, and we used the 
notation: 



0^iE-JL^,y (.M£^ ,21) 



4 The Thirring model with angular momentum. 



Long ago Thirring proposed a very simple model for a star In spite 
of its extreme simplicity, it mimics the main features of gravitating systems 
with surprisingly good success. Thirring considered the model without angular 
momentum. We shall see in this section that, as expected, angular momentum 
drastically changes the behaviour of the system. 
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4.1 The model. 



The model can be described as follows: a set of N particles are confined in 
a spherical volume V. Inside this volume there is a spherical interaction region 
(core) Vb, concentric to V. Particles outside the core ("atmosphere") do not 
interact, and two particles inside the core have a constant attractive potential 
energy. Using the step function 

0Vo(r) = I ifr^Vo 
the interaction energy is given by 



Hr,r') == -__ey„(r)ev„(r'), (23) 

where m is the mass of a particle and G the "gravitational" constant. (We have 
chosen the coupling constant in analogy with a gravitating system). We can 
carry out the N —^ oo limit of the previous section, with Nm = M fixed. The 
potential energy then reads. 



ci>(ri,...,r^) = -^a', (24) 
where a = Jy^ (Pr pir) is the fraction of particles inside Vq. 

4.2 The mean field equations. 

To simplify the computations, we shall consider the model in two dimen- 
sions]^. Notice that in this case the factor 3/2 which appears in front of (3 
in eq. (^) must be substituted by 1. Let us introduce the following nota- 
tion: R and i?o are respectively the radius of the total volume, V ^ and the 
core, Vo, and k = Vo/V^. Using the dimensionless variables e — E/{GM'^) and 
^ = L'^/{2GM^R'^), the mean field equation (|2|) becomes: 

^It should not be very difficult to solve it three dimensions, but it is slightly more cumber- 
some and no qualitative difference is expected. 
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n(r-\ - / exp{fi + Pa + p^ryR^ } if r < i?o 

P^"^' ~ \ exp { /i + /3^rVi?2 } if r > i?o ^ ' 



with 



n a 



K 



e = 2- (26) 



The two self-consistency equations 



a = (frp{r) and 1 - a = / cfrp{r) (27) 

determine a and From them, it is straightforward to derive the foUowing 
equation for a: 



1 _ e-z^e^ 

log a - log(l -a) - Pa ^ -/3^(1 - k) + log ( ^ _ ) • (28) 



Notice that P is positive, by definition. It is possible to eliminate p from (|25D, 
and then the particle distribution reads 



^exp(-/3e^) ifr-G^o 
^(^) = S (i-o)F, ^ / (29) 

V 



exp ifr^^o 



where 



= = -^^^^ . (30) 

Eq. ( p8| ) provides the complete solution of the system. Once solved for a, for 
different values of the energy e and the angular momentum Q, we can compute 
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the entropy S and the temperature T, which, due to the fact that the mass 
distribution maximizes S, is: 



T 



1 



dS[p,e] 
de 



f3. 



(31) 



Notice that a is an order parameter for the coUapsing phase transition, since 
it is the fraction of particles inside the core Vq. Notice also that (p8| ) can have 
more than one solution for some values of e and fi. When this is the case, 
the true solution is that which maximizes the entropy, the others constituting 
metastable states. 



4.3 Limiting cases. 

Let us study the solutions of (^8|) in some limiting cases. 

(a) If ^ 0, then ^ — > and we recover the original Thirring model. There is 
a phase transition separating a low energy collapsing phase, where a ~ 1, 
from a high energy gas phase, with a ^ k. 

(b) If oo, then oo. In this case a ^ for e ^ 451 and a ^ k for 

Afl. No collapse happens. 

(c) For fixed values of and e ^ oo, we have P and then a k. In 
the high energy region, the system behaves as a gas, irrespective of the 
angular momentum, and the mass distribution is homogeneous. 

(d) The most interesting case is the behaviour of the ground state for fixed val- 
ues of il. In the ground state the energy of the system is saturated by the 
rotational energy (which cannot be zero due to the angular momentum) 
and the potential energy, so that no energy remains for thermal motion. 
This is the only microscopical contribution to the thermodynamical state 
when f3 oo, and corresponds to zero temperature. It is easy to see that 
the entropy goes to — oo like — log (3, hence the non-thermal (rotational 
-|- potential) energy must be at its minimum, for, if not, the system can 
change its mass distribution in such a way that some amount of thermal 
energy 1/(3 appears, rising the entropy. In the next subsection we shall 
see what is the structure of the ground state for the different values of fi. 
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4.4 The ground state. 



As we just discussed, when /3 — > oo the non-thermal energy 



~ j^J d-'rr^pir) " T ^^^^ 

must be at its minimum. It is possible to minimize the rotational energy with- 
out modifying the potential energy, just taking a mass distribution with fixed 
a which maximizes the inertial moment. Obviously, this mass distribution cor- 
respond to a fraction a of the particles in the outer layer of the core, and the 
remaining 1 — a fraction in the outer layer of the system. The dimensionless 
inertial moment is then (l/i?^) j <Prr^ p{r) = 1 — a{l — n). Hence, the non- 
thermal energy is 



= l-a{l-n) - T ' 

The minimum of the above function is at either a = 0, a = 1, or at one 
solution of 



— ^ - a = . (34) 



[1 - a{l-K)] 

Eq. (^J) is equivalent to a real polynomial equation of third degree, and there- 
fore has either one or three real solutions. Since the first term in (^^ has a 
pole at a = 1/(1 — k) > 1, there is always a solution at a > 1/(1 — k), which is 
non-physical since a must belong to [0, 1]. Besides this non-physical solution, 
when is small (jsj) has two more solutions. For ^ 1 



a„„, « f7 (1 - k) and a„.„ w ( 1 - y/VL{l - k)\ . (35) 

From the second derivative of (B4) we see that there is an inflection point at 



1 - 

= --. , (36) 

i — K 
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and a„a^ correspond to a local maximum and to a local minimum. We 
always have a^^x < < a„i„. Taking the derivative of ( p4[ ) respect to - 
with a = Q!,„i„ - we easily see that a„i„ is a monotonically decreasing function 
of fl. Analogously, introducing a = a,nin in (p^), taking the derivative respect 
to fl and using (p^), we see that the value of e„t at its local minimum a„i„ is a 
monotonically increasing function of 17. 



Let us study the absolute minimum of (33) for < a < 1. The three 
candidates are a = 0, 1 or a^i„. For small il, from ( |3^ ) we learned that a„i„ > 1, 
and, since e„t(0) = 12 and e„t(l) = U,/ k— 1/2, the absolute minimum correspond 
to a = 1. But when Q. grows e„t(l) increases faster than e„t(0), and will 
reach the physical region, since it decreases. The following three possibilities 
define three critical values of Q.: 



(i) The relative minimum Q!„,i„ reaches the physical region. The critical value 
r^i is given by a^i„{Q.i) = 1. From (34) we obtain fii = k^/(1 — k). 

(ii) The relative minimum a,„i„ ceases to be the absolute minimum. This 
happens for — when e„t(Q!„,i„) — e„t(0). We have the two equations 

f^2(l-K) , „ ^2 aL„ 

Ct,„i„ — -, ; and \l2 ^ 



[l-a„„„(l-K)]2 l-a,„,„ 2 

The solution is ^2 = l/[8 (1 - nf] and = l/[2 (1 - k)]. 

(iii) The value of e„t at a = surpasses that at a = 1. This happens when 
n = rjg, having e„,(0) = e„t(l). Clearly f^g = k/[2(1 - n)]. 



The ground state behaviour depends on the ordering of these critical fi, 
which in its turn depends on k. For k < 1/2 we have ili < ^2,^3- Then 
the ground state is characterized by complete collapse if < fii, with all the 
particle inside the core. For fii < < there is an incomplete, with most - 
but not all - the particle in the core. The value of a varies continuously from 
a = 1 at r^i to a = 1/[2(1 — k)] at ^2. For fl > il2 the ground state has an 
empty core (a — 0). Notice the jump of a at $72- 

For K > 1/2 we have ^3 < i7i and fl^ < ft2- In this case, the ground state 
for f7 < 173 is characterized by a complete collapse, with all the particles inside 
the core. For 17 > 173 no collapse happens (a = 0). In fig. Q we display the 
behaviour of the ground state for both k < 1/2 and k > 1/2. 
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a b 

Figure 1: The structure of the ground state as a function of for a) k < 1/2 
(observe the logarithmic scale in abscissas) and b) k > 1/2. 

4.5 Phase diagram. 

For the discussion of the model at finite temperatures we consider the most 
interesting case, k < 1/2 (atmosphere bigger than the core). The numerical 
solutions displayed in the figures are for k = l/(e'^ — 1) w 0.0524. Exactly as 
for the ground state, there are three different phases, which can be distinguish 
by the behaviour of the order parameter a, which is displayed in fig. |^ as a 
function of the energy, for three values of fi, representative of each phase. Figs. 
^ and ^ display the temperature and the entropy for the same values of Q. It 
can be clearly seen the correlation between the collapsing transition and the 
anomalies in the thermodynamical quantities. Let us discuss in detail what is 
happening on each phase. 

(1) n < fli. There is a complete collapse at low energies, with a ~ 1. The 
gas phase at high energies is separated from the collapsing phase by an 
interval of energies with negative specific heat. Qualitatively, the model 
is similar to the original {fl = 0) Thirring model. The entropy shows the 
characteristic convex intruder in the energy interval where the two phases 
coexist. (Notice that, since the system is not thermodynamically stable, 
van Hove's theorem |l^ does not apply). 

(2) ill < < ri2- At low energies the collapse is not complete, with a < 1. 
For values of near fii the thermodynamical quantities are qualitatively 
equal to those in the low angular momentum phase (negative specific heat). 
When is larger, eq. ( p8| ) has more than one solution. For energies 
larger than a certain value, one of the new solutions becomes the absolute 
maximum of the entropy. This is the origin of the jumps in a and in T that 
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can be seen in figs. |2| and for — 0.05. Notice that the jump occurs 
after a region with negative specific heat. For H. still larger, the jump 
appears before the specific heat becomes negative. However, although in 
the last cases the specific heat is positive everywhere, the entropy still 
has a convex intruder, due to the kink originated by the jump in the 
temperature. 

(3) r2 > O2. There is no collapse at low energies. The specific heat is positive, 
smooth and increases monotonically with the energy. The entropy has no 
convex intruder. 



From the above discussion, it becomes manifest the correlation between the 
collapsing phase transition and the anomalies in the caloric curve (T versus e) 
and in the entropy. The phase diagram in the plane (e, f2 is displayed in fig. 
Ij. Obviously, there is a forbidden region, where the system cannot be since a 
minimal rotational energy is required to keep angular momentum constant. The 
boundary between the forbidden and allowed regions is the T = isotherm. 
Below the forbidden region, there are two phases^ separated by a transition 
region, which is defined through the Maxwell construction. Since botli phases 
are qualitatively different, they must be analytically separated. This means 
that the critical point must be at zero temperature (i.e., on the curve of minimal 
energy). And, in its turn, this implies that the critical point must be at = f^2, 
since this is the point on the zero temperature line where the collapse disappear. 
This is what actually happens, as can be seen in fig. |^. 



For fixed Q, the difference of the energies limiting the transition region de- 
fines the latent heat. We see that the phase transition is first order everywhere, 
with non- vanishing latent heat. From the analysis of eq. (28) when /3 ^ 00 and 
^ ri2 one concludes^ that the latent heat vanishes as O2 — O. Therefore, at 
the critical point, ri2, the latent heat disappears. Since the specific heat is con- 
tinuous at this point, we cannot say that the phase transition becomes second 
order; indeed, the order parameter a jumps at fl2- 



4.6 Discussion. 



Now, we want to address the question of whether the phase diagram de- 
scribed in the previous subsection shares its main properties with those of more 

^The two collapsing phases, Q < f2i and f2i < H < Q2, are distinguished only by the 
behaviour of the order parameter at low energies. The thermodynamical quantities, however, 
behave in a similar way in both cases. Therefore, in the phase diagram we shall not distinguish 
the two low energy phases. 

^ After a careful study of the possible solutions, their entropies, etc., which is not reported 
here in detail. 
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Figure 2: The collapsing order parameter a versus the energy for Q values on 
each of the three phases. To facilitate comparisons, we have redefined the energy 
by a shift, such that the minimum energy is e = for all fl. 



realistic models, or whether these are merely a consequence of the extreme 
simplicity of the model . We can give at least a partial answer. The model 
considered here is unrealistic mainly due to the fact that the particles can only 
collapse in a fixed region. Therefore, there is no place here for more complicated 
regimes like multifragmentation. 

The collapsing phase, at low angular momentum L, with all - or most - the 
particles forming a cluster, is expected to be present in realistic problems. Of 
course, the size of the cluster will increase with angular momentum. For large 
energies both the potential and rotational energy are negligible and we shall 
have a gas phase. There will be then a transition region separating this two 
phases, like in fig. |[ The difference may be, however, that in realistic cases 
these two phases might be connected through intermediate regimes like multi- 
fragmentation. It is clear that for large enough angular momentum the system 
will prefer to collapse into several bodies rather than in a single one. The critical 
point, if any, might be located at non-zero temperature. If this is the case, one 
could expect a true second order phase transition, with divergent specific heat 
and large fluctuations of some local order parameter, like the density. Needless 
to say, these are questions which deserve a careful investigation. 
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Figure 3: The microcanonical temperature versus the energy for three values of 
fl, one on each phase of the system. 



5 Conclusions. 

Let us briefly summarize the results presented in this work. First, we found 
a non-singular expression for the microcanonical distribution of systems with 
conserved, non-vanishing angular momentum. It was obtained by integrating 
out the momenta in the phase space, and it is suitable for Monte Carlo simu- 
lations, as well as for mean field computations. We presented a derivation of 
the mean field equations for systems of classical particles interacting through 
a two body unstable potential, taking into account the conservation of angular 
momentum. 

We discussed, as an application, the properties of the phase diagram of a 
simple model of gravitating particles with non-vanishing angular momentum. 
We found the exact solution in the mean field approach, which served to il- 
lustrate how the phase diagram of a physical system can be modified by the 
conservation of the angular momentum. In the case under interest, the angular 
momentum L ^ ^/fi determines the existence of three phases, separated by two 
critical values Qi and il2. At low energies all particles condense in a small re- 
gion if < 57i For fii < 57 < the collapse at low energy is incomplete with 
most - but not all - particles condensed in a small region. At high energies, 
the system behaves like a gas irrespective of the value of il. The collapsing-gas 
phase transition for < is accompanied by an anomaly in the caloric curve 
(T versus e) , which reflects the fact that the entropy is not concave in an energy 
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Figure 4: The entropy versus the energy for three values of fi, one on each phase 
of the system. The convex intruder for the two smaUer il can be appreciated. 
For = 0.15 the entropy is concave. 



interval. This interval coincides with that where the collapsing-gas transition 
takes place. Through the Maxwell construction, it is possible to define precisely 
this transition region and a latent heat, which shows that the transition can be 
classified as first order. The latent heat disappear continuously at = and 
T = 0. This reflects the fact that both phases, being qualitatively different, 
cannot be analytically connected in the phase diagram. At O > there is nei- 
ther collapse nor anomaly in the thermodynamical quantities. We argued that, 
with more realistic potentials, an intermediate regime between the collapse and 
gas phases - multifragmentation perhaps - should exist. Therefore, it might be 
possible that the system pass smoothly from one phase to the other. Therefore, 
the critical point might be located at non-zero temperature, showing the typical 
behaviour of a second order phase transition. 
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Figure 5: The phase diagram of the model. Notice the logarithmic scale in the 
ordinate axis. The thick line corresponds to the T = isotherm. 
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